Symmetry breaking and manipulation of nonlinear optical modes in an asymmetric 

double-channel waveguide 
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We study light beam propagation in a nonlinear coupler with an asymmetric double-channel 
waveguide, and derive various analytical forms of optical modes. The results show that the 
symmetry-preserving modes in a symmetric double-channel waveguide is deformed due to the asym- 
metry of the two-channel waveguide, and such a coupler supports yet the symmetry-breaking modes. 
The dispersion relations reveal that the system with self-focusing nonlinear response supports the de- 
generate modes, while for self-defocusing medium the degenerate modes do not exist. Furthermore, 
nonlinear manipulation is investigated by launching optical modes supported in double-channel 
waveguide into a nonlinear uniform medium. 
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I. INTRODUCTION 

Propagation of optical waves in waveguide arrays has 
become an important and effective means to investigate 
various optical phenomena which have analogues in many 
fields of physics Special attention has been paid to 
nonlinear surface waves and nonlinear guided waves in 
planar layered structures [3-[l3j and nonlinear couplers 
[Tj-Qil, generation and properties of solitons in nonlin- 
ear waveguide arrays [19M24I ] and other nonlinear periodic 
systems, such as optically induced lattices [25-28]. The 
investigation of light beam propagation in waveguide ar- 
rays attracts increasing attention due to their potential 
applications in all-optical signal processing in fiber optic 
networks and devices, the passive mode-locki ng by u sing 
waveguide array |2{|, and the beam steering [30| - |32| . 

The behavior of light beam propagation in a coupler 
composed by two-channel nonlinear waveguide gained 
particular attention because it can exhibit some of uni- 
versal properties in nonlinear periodic systems and non- 
linear waveguide arrays [33] . It is shown that the coupler 
can support the symmetry-preserving solutions which 
have the linear counterparts and the symmetry-breaking 
solutions without any linear counterparts [3J-[36|], in 
which the spontaneous symmetry-breaking has been ex- 
perimentally demonstrated in optically induced lattices 
with a local double- well potential [34| • 

In this paper, we consider light beam propagation in 
an asymmetric double-channel waveguide with Kerr-type 
nonlinear response, and derive various analytical station- 
ary solutions in detail. It is found that the asymmet- 
ric double-channel waveguide can break the symmetric 
form of the symmetry-preserving modes otherwise in the 
symmetric double-channel waveguide, and such a coupler 
supports the symmetry-breaking modes. We also inves- 



tigate how the type of nonlinear response affects the ex- 
istence and properties of nonlinear optical modes in the 
asymmetric double-channel waveguide. The dispersion 
relation shows that the degenerate modes exist in the 
system with self-focusing nonlinear response, while for 
the coupler with self-defocusing response the degenerate 
modes do not exist. In addition, based on these optical 
modes supported in asymmetric double-channel waveg- 
uide we demonstrate the control and manipulation of op- 
tical modes in different nonlinear media by tuning the 
waveguide parameters. 

The paper is organized as follows. In Section II, the 
model equation describing beam propagation in a double- 
channel waveguide is derived. In Section III, various an- 
alytical forms of optical modes are presented both in self- 
focusing and self-defocusing media. Meanwhile, the dis- 
persion relations between the total energy and the prop- 
agation constant are discussed. In Section IV, we study 
the nonlinear manipulation of optical modes in double- 
channel waveguide. The conclusion is summarized in Sec- 
tion V. 



II. MODEL EQUATION AND REDUCTIONS 

We consider a planar graded-index waveguide with re- 
fractive index 



n(z, x) — F(x) + ri2l(z, x). 



(1) 
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Here the first term on the right hand side presents a two- 
channel waveguide with the different refractive index, 
namely, F(x) — n\\ as — L /2 — Dq < x < —L /2, and 
F(x) = 7ii2 as Lo/2 < x < Lo/2 + Do, otherwise, F(x) = 
no (< nn,ni2), where Do and Lq represent the width 
of waveguide and the separation between waveguides, re- 
spectively; while no, Tin and n\2 denote the refractive 
index of cladding and waveguide, respectively. The sec- 
ond term denotes Kerr-type nonlinearity, I{z, x) is the 
optical intensity, and positive (negative) value of nonlin- 
ear coefficient 712 indicates self-focusing (self-defocusing) 
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FIG. 1: The profile of a two-channel waveguide with different 
refractive index. 



medium. Under slowly varying envelope approximation, 
the nonlinear wave equation governing beam propaga- 
tion in such a waveguide with the refractive index given 
by Eq. ([I} can be written as 



.dtp 1 d 2 i> k [F(x) 
dz 2ko dx 2 



no 



no 



0, (2) 



where ij){z,x) is the envelope function, and fco = 
27mo/A is wave number with A being wavelength of 
the optical source generating the beam. By in- 
troducing the normalized transformations ip(z, x) = 
(fc \n 2 \ L D /n y 1/2 (p(C,0, £ = x/w and C = z/L D with 
Lo = 2fcou>Qj which represents the diffraction length, we 
get the dimensionless form of Eq. ([2]) as follows 



9V 



(3) 



where r\ — n 2 / \n 2 \ — ±1 corresponds to self-focusing (+) 
and self-defocusing (— ) nonlinearity of the waveguides, 
respectively, and V(£) — 2k 2 w 2 [F(wq£) — n ] /n is of 
the form 



v(0 



-L/2 - D < £ < -L/2, 
L/2 < £ < L/2 + D, 
otherwise, 



(4) 



which describes the dimensionless two-channel waveg- 
uide structure with different refractive index, where V\ — 
2k 2 w 2 (nu -n )/n a and V 2 = 2fcoWo("i2 ~n a )/n being 
the modulation depth of the refractive index of the left 
and right waveguide, and L — Lq/wq and D — Dq/wq 
corresponding to scaled separation and width of waveg- 
uide, respectively. Here, we use the typical waveguide 
parameters D = 3.5, L = 5, Vi = 2.525, and vary V\. 
Figure 1 shows the profile of the two-channel waveguide 
structure given by Eq. ©. It should be pointed out that 
such structure can be realized experimentally [38]. It is 
shown that the Eq. ([3]) conserves the total energy flow 

P (0 = l^(C,£)| 2 ^ = Po, where P is the dimen- 
sionless initial total energy. 

Assuming that the stationary solution of Eq. ^ is 
of the form ip(C,£,) = ex P(*/^C): where it(£) is a real 



function, and j3 is the propagation constant, and sub- 
stituting it into Eq. ([3]), we find that the function 
obeys the following nonlinear equation 



(3]u + r/u 3 = 0, 



(5) 



where n = ±1 corresponds to self-focusing (+) and self- 
defocusing (— ) nonlinearity of the waveguides, respec- 
tively. It should be pointed out that Eq. (O not only 
can describe the optical modes in the double-channel 
waveguide structure with the different refractive index, 
but also can describe one-dimensional Bose-Einstein con- 
densate trapped in a finite asymmetry double square 
well potential — V(£). In particular, when Vi = V%, the 
corresponding optical modes in the symmetric double- 
channel waveguide structure have been studied, and the 
results have shown that the coupler not only supports 
symmetry- preserving modes but also symmetry-breaking 
modes 1301. 



III. OPTICAL MODES 

In this section, we will present the analytical solutions 
of Eq. ([S]) with the potential flU for r) = ±1. Generally, 
the solutions of Eq. ([5]) can be constructed in terms of the 
Jacobi elliptic functions depending on the values of the 
variable £, and have the same propagation constants in 
different regions. Within the double-channel waveguides 
of -L/2 - D < £ < -L/2 and L/2 < £ < L/2 + D, the 
solution of Eq. §5§ is the oscillatory function, so it can 
be selected in the form 37] 



Ul (£]A,K,5) =Asn [K£ + 6 



V A 2 \ 
2K 2 



(6) 



with = Vi - K 2 + r)A 2 /2 in the region of —L/2 — 
D < £ < -L/2 and (3 = V 2 - K 2 + r)A 2 /2 in the region 
of L/2 < £ < L/2 + D. In the region of |£| < L/2, 
the solution of Eq. ([5]) has two different Jacobi elliptic 
functions for the symmetric and the antisymmetric case, 
respectively. For the symmetric case, the solution is [37j 



u 2 (£;B,Q,a) = Bnc Q£ + a, I 



T]B 2 

2Q 2 



(7) 



with P = Q 2 + rjB 2 ; and for the antisymmetric case, the 
solution is [37j 



u 2 (£;B,Q,a) = B sc Qf + a, 1 



T]B 2 

2Q2 



(8) 



with ft — Q 2 — T]B 2 /2. It should be noted that those two 
solutions are precise solutions of Eq. (J5J) for one node 
and no node within the region of |£| < L/2. In other 
regions, the solution of Eq. ([5]) is required to tend to 
zero as £ — > ±oo, so it is taken as [36| 
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(9) 
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with /3 > 0. Substituting © into Eq. ©, one obtains 
c = 77/(8/36). 

Note here that although the modulus in the usual Ja- 
cobi elliptic function is restricted from to 1, this prob- 
lem can be solved by the modular transformation such 
that the modulus in the Jacobi elliptic functions given 
by Eqs. can take any positive or negative values in 

our investigations, as shown in Refs. [371. |39( . so those 
solutions do not depend on nonlinearity sign. 

In the following, we show the analytical global solu- 
tions of Eq. ©. With the help of Eqs. ©, © [or ©], 
and © , the solutions of Eq. can be written as 



«(0 



«3(£;&i), 

ui{£;A 1 ,K 1 ,5 1 ), 

U2(£;B,Q,<j), 

ui{^A 2 ,K 2 ,8 2 ), 



£ < -L/2 - D, 
-L/2 -D < £ < -L/2, 
If I < L/2, 

L/2 < £ < L/2 + D, 
£ > L/2 + D. 

( 10 ) 

The continuity conditions for u and du/d^ at the bound- 
aries of £ = ±L/2 and £ = ±(L/2 + D) require 

u 3 {-L/2 - D; bx) = u x {-L/2 - D; A X ,K U 8 X ) , 



du 



-^{-L/2-D-M) = -^(-L/2-D;A 1 ,K 1 ,5 1 ), 



rf£ 



ui {-L/2; A 1 ,K 1 ,S 1 )= u 2 (-L/2; B, Q, a), 



dui 



{-L/2\A U K U S X ) 



du 2 



{-L/2;B,Q,a), 



u 2 {L/2; B, Q, a) = Ul {L/2; A 2 , K 2 , S 2 ), 



du 2 
~d£ 



{L/2;B,Q,a) 



du\ 



{L/2;A 2 ,K 2 ,S 2 ), 



n{L/2 + D; A 2 ,K 2 , S 2 ) = u 3 {L/2 + D; b 2 ), 



du\ 



{L/2 + D; A 2 ,K 2 , S 2 ) = ^{L/2 + D; b 2 ) 



ri£ 



(11) 



with /3 = Vi - K\ + t]A\/2 = V 2 - K'i + i)A\/2, and 
(3 = Q 2 + rjB 2 for the symmetric case given by Eq. © 
or j3 = Q 2 — rjB 2 /2 for the antisymmetric case given by 
Eq. ((HJ. In Eq. (fpQ|) . there are eleven parameters A±, 
K\, 5\, A 2 , K 2 , S 2l B, Q, a, b\, and b 2 , which can be 
calculated by solving numerically Eqs. with the con- 
ditions that the propagation constants in different regions 
should be same. Once those parameters are determined 
numerically, we can obtain the exact optical modes for 
asymmetric double- channel nonlinear waveguides. 

Fig. 2 and Fig. 3 show several different optical 
modes in a nonlinear asymmetric double-channel waveg- 
uide in the self-defocusing medium and the self-focusing 
medium, respectively. These optical modes are induced 
from the symmetry-preserving optical modes in the sym- 
metric double-channel waveguide, where for comparison, 
we also plotted the corresponding symmetry-preserving 
optical modes in the symmetric double-channel waveg- 
uides in the same figure. From Figs. 2 and 3, we found 
that the symmetry of the modes is broken due to asym- 
metry of the two-channel waveguide, and the amplitude 
of the modes in the lower refractive index waveguide is 
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FIG. 2: (Color Online) Various different optical modes ex- 
isting in self-defocusing medium (77 = —1), where the dash- 
dotted red line is the optical mode for Vi = 2.500, the solid 
green line is the optical mode for Vi — 2.525, and the dashed 
blue line is the optical mode for Vi = 2.550. Here /3 = 2.00 
in (a) and (b), /3 = 0.85 in (c) and (d). 
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FIG. 3: (Color Online) Various different optical modes exist- 
ing in self-focusing medium {rj = 1), where the dash-dotted 
red line is the optical mode for Vi = 2.500, the solid green 
line is the optical mode for Vi = 2.525, and the dashed blue 
line is the optical mode for Vi = 2.550. Here /3 = 2.30 in (a) 
and (b), /3 = 1.10 in (c) and (d). 



smaller than that in the higher refractive index waveguide 
for the self-defocusing medium, while for the self-focusing 
medium, the amplitude of the modes in the lower refrac- 
tive index waveguide is larger than that in the higher 
refractive index waveguide. 

We also demonstrate the profiles of optical modes with 
dependence on the propagation constant ft. Fig. 4 and 
Fig. 5 present several corresponding modes shown in 
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FIG. 4: (Color Online) Several optical modes with different 
P in a nonlinear asymmetrical double-channel waveguide for 
the self-defocusing medium (77 = —1). Here the parameters 
are V\ = 2.500 and V 2 = 2.525. 
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FIG. 5: (Color Online) Several optical modes with different 
P in a nonlinear asymmetrical double-channel waveguide for 
the self- focusing medium (r\ — 1). Here the parameters are 
the same as in Fig. 4. 



Figs. 2 and 3 for different propagation constant (3. From 
them, one can see that, for self-defocusing media the pro- 
file of nonlinear modes is shrunk, and the correspond- 
ing amplitude becomes smaller with an increase of the 
propagation constant f3, while for self-focusing case it 
is opposite, namely, the profile of nonlinear modes be- 
comes more prominent and the corresponding amplitude 
becomes larger. This feature can be depicted by the dis- 
persion relations between the total energy Pq and the 
propagation constant /3. As shown in Fig. 8 and Fig. 
9, one can see that the total energy decreases with the 
increase of the propagation constant for self-defocusing 



FIG. 6: (Color Online) Several symmetry-breaking optical 
modes, where the dash-dotted red lines are optical modes for 
V\ = 2.500, the solid green lines are optical modes for Vi = 
2.525, and the dashed blue lines are optical modes for Vi = 
2.550. Here r\ — —1 in (a) and (c), rj = 1 in (b) and (d) and 
P = 1.78 in (a), p = 2.39 in (b), p = 0.66 in (c) and P = 1.15 
in (d). 
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FIG. 7: (Color Online) Several symmetry-breaking optical 
modes with different P in a nonlinear asymmetrical double- 
channel waveguide. Here the parameters are the same as in 
Fig. 6. 



media (see Fig. 8), whereas it is an increasing function 
of the propagation constant for self-focusing case (see Fig. 
9). 

As discussed in Ref. (36|, besides the symmetry- 
preserving optical modes, a double-channel waveguide 
also supports the symmetry-breaking optical modes, and 
the corresponding optical modes in a nonlinear asym- 
metric double-channel waveguide are presented in Fig. 
6, in which we also plotted the corresponding symmetry- 
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FIG. 8: (Color Online) The dependence of the total energy 
Po on the propagation constant /3 for the modes existing in 
self-defocusing medium (77 = —1). Here the parameters are 
Vi = 2.500 and V2 = 2.525. The labeled shadow areas are 
enlarged in corresponding (a)-(b). And the labels 2a, 2b, • ■ • , 
mean the corresponding modes shown in Fig. 2a, 2b, • • • , 
respectively. 
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FIG. 9: (Color Online) The dependence of the total energy Po 
on the propagation constant /3 for the modes existing in self- 
focusing medium (77 = 1). Here the parameters are Vi = 2.500 
and V2 = 2.525. The labeled shadow areas are enlarged in 
corresponding (a)-(d). And the labels 3a, 3b, ■ • • , mean the 
corresponding modes shown in Fig. 3a, 3b, • • • , respectively. 



breaking optical modes in a symmetric double-channel 
waveguide in the same figure for comparison. One can 
find that the optical modes in a nonlinear asymmetric 
double-channel waveguide are almost the same as the 
modes in a symmetric one for the given /3. 

Similarly, the corresponding modes shown in Fig. 6 for 
different propagation constant /3 are presented in Fig. 7. 
It is shown that the profile of nonlinear modes becomes 
shrunk with the increase of the propagation constant (3 
for the self-defocusing case, while for the self-focusing 
case the profile of nonlinear modes becomes more pro- 
nounced. This feature is depicted by the dispersion re- 
lations between the total energy Po and the propagation 
constant j3. It should be pointed out that the modes 
shown in Fig. 6 only exist in a small region, as shown in 
Fig. 8 and Fig. 9. 

From the dispersion relations shown in Figs. 8 and 
9, one can see that for the self-defocusing medium there 
is no intersection between dispersion curves (see Fig. 8 
and the enlarged Figs. 8a and 8b), which indicates that 
there is no degenerate modes existing, and the total en- 
ergy of the mode shown in Fig. 2a is the highest for a 
given propagation constant j3. While for the self-focusing 
medium the dispersion curves can intersect (see the en- 
larged Figs. 9c and 9d), which implies that there exists 
two different modes with the same total energy at the 
intersection point, namely, the degeneracy occurs at the 
intersection point. Note that for our choices of the pa- 



rameters, the intersection points of the dispersion curves 
for the modes shown in Figs. 3a and 3b, and Figs. 3c 
and 3d are about 2.5485 and 1.3595, respectively, and the 
corresponding degenerate modes are shown by the green 
solid curves in Fig. 5. Here, to distinguish the case of 
the intersection, we rotate the dispersion curve for the 
modes shown in Fig. 3a (Fig. 3c) for an angle of 7r/6 
anticlockwise as the center of intersection point and the 
same angle for dispersion curve of the mode shown in 
Fig. 3b (Fig. 3d) but rotate clockwise, as shown in Figs. 
9c and 9d. 



IV. THE NONLINEAR MANIPULATION OF 
OPTICAL MODES IN DOUBLE-CHANNEL 
WAVEGUIDE 

In this section, we will demonstrate the control and 
manipulation of optical modes in reconfigurable non- 
linear media. Thus our interest is to investigate the 
evolution dynamics of optical beams existing in double- 
channel waveguide propagating into a uniform nonlinear 
medium. In this case, the governing equation can be 
generally written as 

M + 1 ^ + M«^o, (12) 
oz Zkq ox* tiq 

where the refractive index change An(z, x) is a function 
of z and x, and An(z, x) = n{z, x) — tiq. Here we assume 
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FIG. 10: (Color Online) The evolution of optical modes shown 
in Fig. 2 into the self-focusing Kerr medium without channels, 
where rj = 10 in (a) and (b), and rj = 20 in (c) and (d). Here 
the parameters are Vi = V% = 2.525 and £o = 10. The labels 
(a), (b), (c), and (d) mean the corresponding modes shown in 
Fig. 2a, 2b, 2c, and 2d, respectively. 



FIG. 11: (Color Online) The evolution of optical modes shown 
in Fig. 3 into the self-focusing Kerr medium without channels, 
where rj = 5 in (a) and (b), and rj — 10 in (c) and (d). Here 
the parameters are V% = V% = 2.525 and fo = 10. The labels 
(a), (b), (c), and (d) mean the corresponding modes shown in 
Fig. 3a, 3b, 3c, and 3d, respectively. 



that when < z < Zo, n(z,x) is in the form of Eq. 
([T]), An(z, x) = F{x) + n,2l{z, x) — no; while for z > Zq, 
An(z, x) = n' 2 I{z, x) — tiq. Here, ri2 and n' 2 are the Kerr 
nonlinear coefficients of different media in the region of 
< z < Zq and z > Zq, respectively. Thus, when < 
z < Zo, namely in the region of < ( < Co, Eq. (O 
can be normalized to Eq. ([3]), where £ = z/L]j and 
Co = Zq/Lu with Ld = 2fcoWg- At the same time, in the 
region of z > Z , namely C > Co, Eq. (fT2")l is reduced to 
the dimensionless form as follows 



■ di P , d2( f , Mi 2 n 



(13) 



where rj' = n' 2 / \n 2 \- Note that Eq. p^|) is different from 
Eq. in which Eq. includes a potential function 
V(£) given by Eq. (gj, while Eq. (IT3l) does not include 
and can describe the dynamics of beams in Kerr media 
without any refractive index modulation. 

In the following analysis, optical beams of different 
modes existing in double-channel waveguide are injected 
into the uniform nonlinear medium after propagating 
diffraction length of Co in double-channel waveguide. 

First we consider the situation that optical beams 
are injected from symmetrical double-channel waveguide. 
For the optical modes shown in Fig. 2, which exist in 
self-dcfocusing medium (r) = —1), the numerical sim- 
ulations show that when -q' < 0, these optical modes 
arc diffracted quickly after entering into a uniform Kerr 
medium. However, when rj' > and is large enough, the 
evolution of optical beams exhibits different scenarios in 
the Kerr medium without any channels, as shown in Fig. 



10. Similarly, for the optical modes shown in Fig. 3, 
which exist in self- focusing media (rj = 1) with double- 
channel waveguide, as shown in Fig. 11, our numerical 
simulations show that the evolution of optical modes ex- 
hibit almost similar properties with that in Fig. 10. 

From Figs. 10 and 11, one can see that when the 
optical modes existing both in self-defocusing and self- 
focusing media with double- channel waveguide are in- 
jected into the self-defocusing medium without any chan- 
nels, the beams should be diffracted quickly. However, 
when the optical modes are injected into the self-focusing 
medium without any channels and the corresponding 
nonlinear coefficient r/ is enough larger, the beams could 
be manipulated effectively. In this situation, when the 
optical modes shown in Fig. 2a and Fig. 3a are in- 
jected into self-focusing medium without any channels 
the beams appear to attract and repel each other peri- 
odically, as shown in Figs. 10a and 11a. While when 
the modes shown in Figs. 2b-2d and Figs. 3b-3d are 
injected into self-focusing medium without any channels, 
the beams repel each other, as shown in Figs. lOb-lOd 
and Figs, llb-lld. Note that the escape angle are the 
same for the beams in Fig. 2b (Fig. 3b) due to the 
symmetry of the double-channel waveguide. 

The evolution of optical modes shown in Fig. 6 is pre- 
sented in Fig. 12, in which Figs. 12a and 12c (Figs. 
12b and 12d) demonstrate the evolution of optical modes 
in the self-focusing medium without any channels with 
initial input beams injected from self-defocusing (self- 
focusing) medium with double-channel waveguide. As 
shown in Fig. 12a and 12b, one can see that optical 
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FIG. 12: (Color Online) The evolution of optical modes shown 
in Fig. 6 into the self-focusing Kerr medium without chan- 
nels, where rj — 10. Here the parameters are Vi = Vi = 2.525 
and Co = 10. The labels (a), (b), (c), and (d) mean the corre- 
sponding modes shown in Fig. 6a, 6b, 6c, and 6d, respectively. 



It should be pointed out that these results only take 
into account the optical modes existing in the symmetri- 
cal double-channel waveguide. Then, one will naturally 
ask what is the influence of the asymmetrical double- 
channel on the evolution of optical modes. In order to 
understand this question, we launch optical beams from 
an asymmetrical double-channel waveguide into the self- 
focusing medium to observe their evolution by tuning the 
depth of left channel of the waveguide for a fixed depth 
of the right channel, namely, the value of V\ varies from 
2.500 to 2.550 for Vi — 2.525. As an example, we demon- 
strated the evolution dynamics for the modes shown in 
Fig. 2b and Fig. 3b. In Fig. 13, we present the de- 
pendence of the energy sharing, the ratio of the energy 
carried by each component in the mode over the total en- 
ergy, and the escape angle, the angle of the each peak in 
the mode and the propagation direction on the value 
of V\. As shown in Fig. 13, one can see that the energy 
carried by each beam is different due to the asymmetry 
of the double-channel waveguide [shown in Figs. 13a and 
13b]. Also, one can see clearly that the escape angles of 
the two beams take different values with the change of the 
value Vi, which means that the beams can be controlled 
by tuning the depth of the left channel of the waveguide. 




t , 1 , 1 , 1 , 1 , 1 0.10 I ' 1 ' 1 ' 1 ' 1 ' — I 

250 2.51 2.52 2.53 2.54 2 55 2.50 2.51 2.52 2.53 2.54 2.55 



FIG. 13: (Color Online) Energy sharing (row 1) and escape 
angles (row 2) of optical beams as a function of parameter 
Vi . In all cases, the solid-blue and the dashed-red curves cor- 
respond to the left and right beams, respectively. Here the 
parameter are Vi = 2.525, rj — 10 in (a) and (c) correspond- 
ing to the modes shown in Fig. 2b (namely Fig. lib) and 
rj' — 5 in (b) and (d) corresponding to the modes shown in 
Fig. 3b (namely Fig. 12b). 

beams with a single hump can be compressed effectively, 
and as shown in Figs. 12c and 12d the optical modes 
with two peaks are separated during the evolution due 
to the repulsive interaction force resulted from the phase 
difference between the two peaks. 



V. CONCLUSION 

We have studied light beam propagation in an asym- 
metric double-channel waveguide in the form of a nonlin- 
ear coupler. A family of analytical solutions with sym- 
metric and antisymmetric forms has been obtained for 
both self-focusing and self-defocusing nonlinear media, 
and the dispersion relations between the total energy 
and the propagation constant has been discussed in de- 
tail. Our results reveal that the system with self-focusing 
nonlinear response supports the degenerate modes, while 
for self-defocusing medium the degenerate modes do not 
exist. In addition, we explored new ways to steer op- 
tical modes propagating from double-channel waveguide 
into a uniform self-focusing medium. The compression 
of beam with single hump and split of beams with two 
humps have been demonstrated by tuning the depth of 
the channel of the waveguide. These properties may be 
applied in practical optical devices, and be useful for op- 
tical processing, optical switching or optical routing. 

VI. ACKNOWLEDGEMENT 

The authors acknowledge fruitful discussions with Pro- 
fessor Yuri Kivshar. This research is supported by the 
National Natural Science Foundation of China grant 
61078079, the Shanxi Scholarship Council of China grant 
2011-010, and the Australian Research Council. 



8 



[1] S. Longhi, Laser Photon. Rev. 3, 243 (2009). 

[2] W. J. Tomlinson, Opt. Lett. 5, 323 (1980). 

[3] V. K. Fedyanin and D. Mihalache, Z. Phys. B47, 167 

(1982); D. Mihalache, R. G. Nazmitdinov, and V. K. 

Fedyanin, Physica Scripta 29, 269 (1984); D. Mihalache, 

D. Mazilu, and H. Totia, Physica Scripta 30, 335 (1984). 
[4] N. N. Akhmediev, Zh. Eksp. Teor. Fiz. 83, 545 (1982) 

[Sov. Phys. JETP 56, 299 (1982)]. 
[5] G. I. Stegeman, C. T. Seaton, J. Chilwell, and S. D. 

Smith, Appl. Phys. Lett. 44, 830 (1984). 
[6] N. N. Akhmediev, V. I. Korneyev, and Y. V. Kuzmenko, 

Zh. Eksp. Teor. Fiz. 88, 107 (1985). 
[7] F. Lederer and D. Mihalache, Solid State Commun. 59, 

151 (1986); D. Mihalache, D. Mazilu, and F. Lederer, 

Opt. Commun. 59, 391 (1986). 
[8] D. Mihalache, G. I. Stegeman, C. T. Seaton, E. M. 

Wright, R. Zanoni, A. D. Boardman, and T. Twardowski, 

Opt. Lett. 12, 187 (1987). 
[9] D. Mihalache, M. Bertolotti, and C. Sibilia, Prog. Opt. 

27, 229 (1989). 

[10] U. Langbein, F. Lederer, T. Peschel, U. Trutschel, and 
D. Mihalache, Phys. Reports 194, 325 (1990). 

[11] A. D. Boardman, P. Egan, U. Langbein, F. Lederer, and 
D. Mihalache, in "Nonlinear surface electromagnetic phe- 
nomena", Edited by H. E. Ponath and G. I. Stegeman, 
North-Holland, Amsterdam, 1991, pp. 73-287. 

[12] F. Lederer, L. Leine, R. Muschall, T. Peschel, C. 
Schmidt-Hattenberger, T. Trutschel, A. D. Boardman, 
C. Wachter, Opt. Commun. 99, 95 (1993). 

[13] A. R. Davoyan, I. V. Shadrivov, and Y. S. Kivshar, Opt. 
Express 16, 21209 (2008). 

[14] D. Marcuse, Journal of Lightwave Technology 5, 113 
(1987). 

[15] C. C. Yang, Opt. Lett. 16, 1641 (1991). 

[16] C. C. Yang and A. J. S. Wang, IEEE Journal of Quantum 

Electronics 28, 479 (1992). 
[17] Y. Chen, Optical and Quantum Electronics 24, 539 

(1992). 

[18] K. Yasumoto, H. Maeda, and N. Maekawa, Journal of 

Lightwave Technology 14, 628 (1996). 
[19] D. N. Christodoulides, F. Lederer, and Y. Silberberg, 

Nature (London) 424, 817 (2003). 
[20] F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. As- 

santo, M. Segev, Y. Silberberg, Phys. Rep. 463, 1 (2008). 
[21] Y. V. Kartashov, V. A. Vysloukh, and L. Torner, Prog. 



Opt. 52, 63 (2009). 
[22] H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. 
Boyd, and J. S. Aitchison, Phys. Rev. Lett. 81, 3383 
(1998); D. Mandelik, H. S. Eisenberg, Y. Silberberg, R. 
Morandotti, and J. S. Aitchison, Phys. Rev. Lett. 90, 
053902 (2003). 

[23] M. J. Ablowitz and Z. H. Musslimani, Phys. Rev. Lett. 

87, 254102 (2001). 
[24] Y. Lahini, E. Frumker, Y. Silberberg, S. Droulias, K. 

Hizanidis, R. Morandotti, and D. N. Christodoulides, 

Phys. Rev. Lett. 98, 023901 (2007). 
[25] J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. 

Christodoulides, Nature (London) 422, 147 (2003); J. W. 

Fleischer, T. Carmon, M. Segev, N. K. Efremidis, and D. 

N. Christodoulides, Phys. Rev. Lett. 90, 023902 (2003). 
[26] D. Neshev, E. Ostrovskaya, Y. Kivshar, and W. Kro- 

likowski, Opt. Lett. 28, 710 (2003). 
[27] Y. V. Kartashov, V. A. Vysloukh, and L. Torner, Phys. 

Rev. Lett. 93, 153903 (2004). 
[28] Y. P. Zhang and B. Wu, Phys. Rev. Lett. 102, 093905 

(2009). 

[29] J. Proctor and J. Kutz, Opt. Lett. 30, 2013 (2005). 
[30] A. B. Aceves, T. Peschel, R. Muschall, F. Lederer, S. 

Trillo, and S. Wabnitz, Phys. Rev. E 53, 1172 (1996). 
[31] H. F. Zhang, J. Jia, S. T. Jia, and L. Li, Opt. Commun. 

281, 4130 (2008). 
[32] Z. Xu, Mario I. Molina, and Yuri S. Kivshar, Phys. Rev. 

A 80, 013817 (2009). 
[33] Y. V. Kartashov, B. A. Malomed, and L. Torner, Rev. 

Mod. Phys. 83, 247 (2011). 
[34] P. G. Kevrekidis, Zhigang Chen, B. A. Malomed, D. J. 

Frantzeskakis, and M. I. Weinstein, Phys. Lett. A 340, 

275 (2005). 

[35] Ze'ev Birnbaum and B. A. Malomed, Physica D 237, 
3252 (2008). 

[36] J. F. Jia, Y. P. Zhang, W. D. Li, and L. Li, Opt. Com- 
mun. 283, 132 (2010). 

[37] W. D. Li, Phys. Rev. A 74, 063612 (2006). 

[38] T. Pertsch, P. Dannberg, W. Elnein, A. Brauer, and F. 
Lederer, Phys. Rev. Lett. 83, 4752 (1999). 

[39] P. F. Byrd and M. D. Friedman, Handbook of Elliptic 
Integrals for Engineers and Scientists, 2nd ed. Springer- 
Verlag, New York, 1971. 



